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Application of the heavy quark expansion: \Vub\ and spectral moments 

M Luke 

Department of Physics, University of Toronto 

The use of inclusive B decays to determine \Vub\, IK-sl, and heavy quark matrix elements via spectral moments is discussed. 



1 \Vut\ 

A precise and model independent determination of the 
magnitudes of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements Vut is important for testing the Standard 
Model at B factories via the comparison of the angles and 
the sides of the unitarity triangle. 

\V„b\ is notoriously hard to measure model-independently. 
The first extraction of \Vub\ from experimental data re- 
lied on a study of the lepton energy spectrum in inclusive 
charmless semileptonic B decay [[0, a region in which 
(as will be discussed) the rate is highly model-dependent. 
\Vub\ has also been measured from exclusive semileptonic 
B ptv and B — > jitv decay [ 2 1. These exclusive deter- 
minations also suffer from model dependence, as they rely 
on form factor models (such as light-cone sum rules [ 3|) 
or quenched lattice calculations at the present time (for a 
review of recent lattice results, see [RH). 

In contrast, inclusive decays are quite simple theoretically, 
and if it were not for the huge background from decays to 
charm, it would be straightforward to determine \ V„b\ from 
inclusive semileptonic decays. Inclusive B decay rates can 
be computed model independently in a series in Aqcd/'«/j 
and as(mb) using the heavy quark expansion (HQE) [ 
E |6l El m. At leading order, the B meson decay rate is 
equal to the b quark decay rate. The leading nonpertur- 
bative corrections of order A^^^jm^^ are characterized by 
two heavy quark effective theory (HQET) matrix elements, 
usually called A\ and /I2. These matrix elements also oc- 
cur in the expansion of the B and B* masses in powers of 
Aqcd/^w^, 



niB(B*) - mb+ A- 



M+3(-l)A2 
2mh 



+ . . 



(1) 



Similar formulae hold for the D and D* masses. The pa- 
rameters A and Ai are independent of the heavy b quark 
mass, while there is a weak logarithmic scale dependence 
in /I2. The measured B* - B mass splitting fixes A2(mb) = 
0.12 GeV^, while (or, more precisely, a well-defined 
short-distance mass such as [ 9. .10.1 ) and Ai may be 
determined from other physical quantities (as will be dis- 
cussed in the second part of this talk). Since the parton 



level decay rate is proportional to m^, the uncertainty in 
nib is a dominant source of uncertainty in the relation be- 
tween B — » X„{v and \Vub\', an uncertainty in nib of 50 MeV 
corresponds to a ~ 5% determination of IV„hirl^ fm . 

Unfortunately, the B — > Xu(v rate can only be measured 
imposing severe cuts on the phase space to eliminate the 
-100 times larger B —> X^tv background. Since the pre- 
dictions of the OPE are only model independent for suf- 
ficiently inclusive observables, these cuts can destroy the 
convergence of the expansion. This is the case for two 
kinematic regions for which the charm background is ab- 
sent and which have received much attention: the large lep- 
ton energy region. Eg > (m^ - m^)/2mB, and the small 
hadronic invariant mass region, nix < nio r il2lll3|[T^ll5l . 

The poor behaviour of the OPE for these quantities is 
slightly subtle, because in both cases there is sufficient 
phase space for many different resonances to be produced 
in the final state, so an inclusive description of the decays 
is still appropriate. However, in both of these regions of 
phase space the B —> Xu(v decay products are dominated 
by high energy, low invariant mass hadronic states. 



nib, ni^ -nijj- Aqcd'W/p '«i 



(2) 



In this region the differential rate is very sensitive to the 
details of the wave function of the b quark in the B meson 
r il6l . This can be seen simply from the kinematics. A b 
quark in a B meson has momentum 



(3) 



where v'^ is the four-velocity of the quark, and is a 
small residual momentum of order Aqcd- If the momen- 
tum transfer to the final state leptons is q, the invariant mass 
of the final state hadrons is 



— inihV + k — q) = {ifibV - q) 
+ 2k ■ {nibV -q) + O(Aqcd)- 



(4) 



Over most of phase space, the second term is suppressed 
relative to the first by one power of Aqcd /'"A, and so may 
be treated as a perturbation. This corresponds to the usual 
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OPE. However, in the region (|2j Ex is large and nix is 
small, rnhvi^ -q'' - (Ex, 0,0, Ex) + 0( Aqcd) is almost light- 
like, and the first two terms are the same order. 



mj( - {nihV — q)^ + 2Exk+ + . . . , k+ = ko + k-},. 



(5) 



The differential rate in this "shape function region" is 
therefore sensitive at leading order to the wave function 
f{k+) which describes the distribution of the light-cone 
component of the residual momentum of the b quark. f{k+) 
is a nonperturbative function and cannot be calculated an- 
alytically, so the rate in the region (|2} is model-dependent 
even at leading order in Aqcd/'m/)- 

Near the endpoint, the OPE for the Ei spectrum has the 
form (where y - 2E(/mb) 



dr 

dy 



29(1 -y)- -^S'(l -y)- ^6"(l -y) + 
3ml 9ml 

- ^6(l-y)-^6(l-y) + ... 



3m 



m^ 



(6) 



For I - y ~ Aqcd /mb, the terms on the first line are all 
parametrically 0(1), and so must be summed to all orders, 
and the result may be written as a convolution of f(k+) with 
the parton-level rate [ 16|. The terms of the second line are 
less singular near y - I, and so correspond to subleading 
effects, which will be discussed later. 

The situation is illustrated in Fig. H^a-b), where the lep- 
ton energy and hadronic invariant mass spectra are plotted 
in the parton model (dashed curves) and incorporating a 
simple one-parameter model for the distribution function 
(solid curves) l .lli 



32 



where 



.;c)2e-;?(i--') 6>(1 -x) 



(7) 



x=— , A = 0.48GeV. 
A 



(8) 



The differences between the curves in the regions of inter- 
est indicate the sensitivity of the spectrum to the precise 
form of /(fc+). In both curves, the unshaded side of the ver- 
tical line denotes the region free from charm background. 
Because ~ Aqcd^mbj the integrated rate in this region is 
very sensitive to the form of f(k+), complicating the issue 
of determining \ V,ib\ model-independently. 



1.1 Optimized Cuts 

One solution to the problem of sensitivity to nonperturba- 
tive effects is to find a set of cuts which eliminate the charm 
background but do not destroy the convergence of the OPE, 
so that the distribution function f(k+) is not required. In 
Ref. [ 181 it was pointed out that this is the situation for a 
cut on the dilepton invariant mass. Decays with 



> (ms - mof 



(9) 



must arise from b ^ u transition. Such a cut forbids the 
hadronic final state from moving fast in the B rest frame, 
and simultaneously imposes mx < mo and Ex < mo- Thus, 
the light-cone expansion which gives rise to the shape func- 
tion is not relevant in this region of phase space [ [2] 1191 . 
The effect of convoluting the q^ spectrum with the model 
distribution function in Eq. Q is illustrated in Fig. Qc). 
The region selected by a q^ cut is entirely contained within 
the cut, but because the dangerous region of high en- 
ergy, low invariant mass final states is not included, the 
OPE does not break down. The price to be paid is that the 
relative size of the unknown Ag^p/m^ terms in the OPE 
grows as the q^ cut is raised. Equivalently, as was stressed 
in r 1201 . the effective expansion parameter for integrated 
rate inside the region (|9j is Aqcd/'Wc, not Aqcd/w^- In ad- 
dition, the integrated cut rate is very sensitive to mi,, with 
a iSOMeV error in mt, corresponding to a ~ +10% uncer- 
tainty in |y„fc| [HOI ED- 

A further important source of uncertainty arises from weak 
annihilation (WA) graphs [ 22|. WA arises at (9(AQ^jj/m^) 
in the OPE, but is enhanced by a factor of ~ 167r^ because 
there are only two particles in the final state compared with 
b —> utvc- Because WA contributes only at the endpoint of 
the q^ spectrum, it is independent of and mcut: 



dr„ 



d^2 



(B2 - Bx)5(q^ - ml). 



(10) 



B\ and B2 are matrix elements which are equal for both 
charged and neutral B's under the factorization hypoth- 
esis, and so the size of the WA effect depends on the 
size of factorization violation. Assuming factorization 
is violated at the 10% level gives a corresponding uncer- 
tainty in \Vub\ from a pure q^ cut of ~ 10% [ 22|; how- 
ever, this estimate is highly uncertain, being proportional 
to \6n^ X (factorization violation). In addition, since the 
contribution is fixed at maximal q^, the corresponding un- 
certainty grows as the cuts are tightened, reducing the inte- 
grated rate. 

These uncertainties may be reduced by considering more 
complicated kinematic cuts: in [ 21 1 it was proposed that 
by combining cuts on both the leptonic and hadronic in- 
variant masses the theoretical uncertainty on |y„;,| could be 
minimized. For a fixed cut on mx, lowering the bound on 
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Figure 1. The shapes of the lepton energy, hadronic invariant mass and leptonic invariant mass spectra. The dashed curves are the b 
quark decay results to OCff.,), while the solid curves are obtained by convoluting the parton-level rate with the model distribution function 
in Eq. 0. The unshaded side of the vertical lines indicate the region free from charm background. 
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Figure 2. The effect of the model structure function on 
G(^^„,,m(.ut) as a function of q^^^ for n\at = 1.86 GeV (solid line), 
1.7 GeV (short dashed line) and 1.5 GeV (long dashed line). 



and subleading twist corrections, this spectrum is directly 
proportional to the structure function. 



dT 

dEy 



Gl\V>tV:/a\q 



eff,2^5 



327r4 



f(Ey). 



(12) 



Thus, combining data on B ^ Xsj with data from B 
Xutv, one can eliminate the dependence on the structure 
function and therefore determine \Vub\ with no model de- 
pendence at leading order [tl6i|24|. At tree level, the rela- 
tion is 



(13) 



increases the cut rate and decreases the relative size of 
the l/'w^ terms (including the WA terms), while only in- 
troducing a small dependence on f{k+). Since this depen- 
dence is so weak, a crude measurement of f{k+) suffices 
to keep the corresponding theoretical error negligible. The 
sensitivity to mt, is also reduced. 

Defining the function G(q^^^^, /Wcut) by 

j^^^^3 G(^,„„mcut), (11) 

the dependence of G{q\^^^, nicut) on fik+) for various cuts is 
illustrated in Fig.|2 The estimated uncertainty from other 
sources is given for a variety of cuts in Tabled Since the 
uncertainty in \Vub\ is half that of G(ql^^, HJcut), we see that, 
depending on the cuts, theoretical errors at the 5-10% 
level are possible. 

1.2 /(fc+) and subleading corrections 

Alternatively, one can reduce the theoretical uncertainty in 
\Vub \ by measuring the universal structure function fik+) in 
some other process f 16 .23J . The best way to measure the 
structure function f(k+) is from the photon energy spec- 
trum of the inclusive decay B — > Xsj. Up to perturbative 



where 

r„{Ec) = \ dEt-^ (14) 

Je, dEt 

2 <^s 
Y,{Ec) = — dEy{Ey - Ec)—— 

nib Je, dEy 

and 5{Ec) contains terms suppressed by OiKQcoImb)- An 
analogous relation holds for the hadronic invariant mass 
spectrum [El 1141 . In addition to higher twist effects, 
there are perturbative corrections to jl3> . Most important 
of these are the parametrically large Sudakov logarithms, 
which have been summed to subleading order 1 1241 . In ad- 
dition, contributions from additional operators which con- 
tribute to B — > Xsj have been calculated [ 25J. The CLEO 
collaboration [ 26 1 recently used a variation of this ap- 
proach to determine \Vub\ from their measurements of the 
B — > Xsj photon spectrum and the charged lepton spec- 
trum in B — > Xuiv(. 

The subleading corrections to ( I13t contained in 6(Ec) have 
only recently been studied [ 27, ^8* "29*, '30]. These are 
analogous to higher twist effects in DIS, and there are two 
separate effects, each of which is large. 

First of all, the 0(1 /nib) corrections happen to have a large 
numerical prefactor. This is easiest to see by looking at the 
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Table 1. G(ql^^f,m^at), as defined in Eq. <1 1> . for several different choices of (^^m, mem), along with the uncertainties. The fraction of 
B — » Xiifv events included by the cuts is 1.21 G(^^„,, OTcui). AstructG gives the fractional effect of the structure function f(k+) in the simple 
model 0; we do not include an uncertainty on this in our error estimate. The overall uncertainty AG is obtained by combining the other 
uncertainties in quadrature. The two values correspond to Amj^* = ±80MeV and ±30MeV. The uncertainty in \ V„i,\ is half of AG. 



Cuts on (q^, nix) 


Giql^v '"cut) 


■^structG 


ApeitG 


+80/30 MeV 




AG 


Combined cuts 


6GeV^1.86GeV 


0.38 


-4% 


4% 


13%/5% 


6% 


15%/9% 


SGeV^, 1.7 GeV 


0.27 


-6% 


6% 


15%/6% 


8% 


18%/12% 


llGeV^1.5GeV 


0.15 


-7% 


13% 


18%/7% 


16% 


27%/22% 


Pure cuts 


(niB - nioY, niD 


0.14 




15% 


19%/7% 


18% 


30%/24% 


(niB - mo-)^,mD- 


0.17 




13% 


17%/7% 


14% 


26%/20% 



OPE for the lepton spectrum in semileptonic b ^ u decay, 
Eqn. (|6j. The terms in the first line are universal, and so 
are the same for B Xsj and B — > X„{vc decays. The 
subleading terms on the second line are not universal, and 
sum to subleading distribution functions [ 27|. Note that 
the subleading A2 term has a coefficient of 11 , whereas the 
corresponding coefficient in the B Xsj is 3, giving an 
0{AQCD/iiib) correction which is enhanced by a factor of 
8 over the naive dimensional estimate. Since this is just 
the first term of an infinite series, one cannot immediately 
determine the size of the subleading correction, but for a 
simple model the corresponding shift of \Vub\ is plotted in 
Fig. |3] For a charged lepton cut of 2.3 GeV, this corre- 
sponds to a ~ 15% shift in the extracted value of \Vub\- 
While this is a substantial shift, it was argued in [ 30 1 that 
the magnitude of this shift is quite insensitive to the model 
chosen for the subleading distribution function, and so the 
coiTesponding uncertainty in \Vub\ is much smaller than the 
overall shift of ~ 15%. 

A second source of uncertainty arises because of the WA 
graphs discussed in the previous section. In the region 
near y = 1, the WA graph is the first term of an infinite 
series which resums into a sub-subleading (relative order 
(1/m^)) distribution function [ 29|. As before, the size of 
the WA contribution is difficult to determine reliably; the 
authors of [ i29i estimate the corresponding uncertainty in 
\Vub\ to be at the ~ 10% level (with unknown sign) for a cut 
E( > 2.3 GeV. For both subleading effects, the fractional 
uncertainty in \Vub \ is reduced considerably as the cut on Ef 
is lowered below 2.3 GeV. 

1.3 Summary for \Vub\ 

We are left in the fortunate situation of having a number 
of theoretically clean methods of determining \Vub\ from 
inclusive decays, each of which has advantages and disad- 
vantages which are summarized in Tabled Because the 
different techniques have different sources of uncertainty, 
agreement between these methods (combined with future 
unquenched lattice predictions for B n(v decays) will 




L ^ ^ J J 

2 2.1 2.2 2.3 2.4 
Ec(GeV) 

Figure 3. The subleading twist corrections to the \V„b\ relation 
<13> as a function on the lepton energy cut, using the simple 
model described in the text, from Ref. [ 1281 . The vertical line 
corresponds to the kinematic endpoint of the semileptonic b ^ c 
decay. 



give evidence that the different sources of uncertainty have 
been correctly estimated. 

In the future, experimental measurements can help reduce 
the theoretical errors in a number of ways: 

• better determinations of nib (through, for example, 
moments of B decay distributions) can reduce the 
largest single source of uncertainty for determina- 
tions using optimized cuts, 

• the size of WA effects may be tested by compar- 
ing Z)° and Ds semileptonic decays, or by extracting 
\Vub\ from B'^ and B^ decays separately, 

• improved measurements of the B Xsj spectrum 
will give an improved determination of fik+), and 

• studying the dependence of the extracted value of 
\Vub\ as a function of the lepton cut Ee can test the 
size of the subleading twist terms in il3\ . 
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Table 2. Comparison of different kinematic cuts for the determination of | V„6| from inclusive decays. 



Cut 


% of rate 


Good 


Bad 








- depends on f{k^) (and subleading 


tf > — s — - 




- simplest to measure 


corrections) 

- WA corrections may be substantial 

- reduced phase space - duality issues? 


nix < nin 


~ 70% 


- lots of rate 


- depends on j(k ) (and subleading 
corrections) 








- very sensitive to mt, 


> (niB - niof- 


-20% 


- insensitive to f{k^) 


- WA corrections may be substantial 

- effective expansion parameter is 

Aqcd/wJc- 






- insensitive to f{k^) 




"Optimized cuts" 


up to 


- lots of rate 


- less rate than pure nix cut. 




-45% 


- can move cuts away from 
kinematic limits and still have 
small uncertainties 


and more complicated to measure 



2 Spectral Moments 

Since differential rates may be computed in the HQE as a 
power series in Q'.s(mi) and Aqcd/w/;, an unlimited number 
of spectral moments may be computed. Different moments 
have different dependence on the nonperturbative parame- 
ters of the HQE, so a simultaneous fit to multiple moments 
allows these parameters to be determined experimentally. 
In addition, consistency between different observables pro- 
vides a powerful check of the validity of the HQE to in- 
clusive decays. Moments which have received particular 
attention are moments of the charged lepton energy spec- 
trum [ 3 1 1 and hadronic invariant mass spectrum [ 32 1 in 
B — > Xciv decays, and moments of the photon energy spec- 
trum in B Xsj decays [ 33 1. 

By comparing the first moment of the photon spectrum in 
B — > Xsy with the first moment of the hadronic invari- 
ant mass spectrum, the CLEO collaboration [ 1341 deter- 
mined A = 0.35 + 0.07 + 0.10 GeV and = -0.236 + 
0.071 + 0.078 GeV^ (where the first error is experimental 
and the second theoretical). More recently, additional mo- 
ments have been measured by CLEO['35'|, BABAR [ |3^ 
and DELPHI [ 37|, providing enough constraints to per- 
form a global fit to the HQE including terms of order 1 /m^. 
Two such global fits have been recently performed. Ref. [ 
|38J found 



AS 



V, 



cb 



4.75 + 0.10 GeV 
(40.8 + 0.9) X 10 



(15) 



-3 



while Ref. [ 1391 found, from just the DELPHI moments. 



nibil GeV) = 4.59 + 0.08 + 0.01 GeV 
m,(lGeV) = 1.13 + 0.13 +0.03 GeV 
Vch = (41.1 + 1.1) X 10"^ 



(16) 



(coiTesponding to m^'^ = 4.69 GeV.) Ref. [ |39 1 also per- 
formed a fit using the pole mass scheme. Matrix elements 
arising up to (9(1 /m^) in the HQE were also determined 
from each of the fits. 

In the approach of Ref. [[38j, the charm quark mass is de- 
termined by the heavy quark relation 

/I 1 \ 

nib-nic - niB -nio- Ai \ — (17) 



+(Pi -Ti -T3) 



Inic Inij, 

Cj i_ 

\ml Ami, 



+ 



(where niBip) = (mB(D) + ^niB'{D'))l^ is the spin-averaged 
meson mass, and p\ and the t, 's are matrix elements of or- 
der Aqpjj). In the approach of Ref. [ 39 1, the charm quark 
is not treated as heavy, and its mass is taken to be a free pa- 
rameter, to be fit from the moments. The second approach 
has the advantage that it does not require an expansion in 
Aqcd/wc- It however has the disadvantage that if the c 
quark is not treated as heavy, a systematic determination of 
the parameters A2 and p2 (arising at 0(1 /m-) and 0(1 /m^), 
respectively) from the D-D* and B - B* mass splittings 
can no longer be performed. 

It is amusing to note that setting experimental errors in the 
global fit of Ref. [ 38 1 to zero, one obtains the uncertainties 
6(\Vcb\ = ±0.35 X lO-\ 6(mb) = +35 MeV, which give 
an indication of the limiting theoretical uncertainty in the 
analysis. 

One can also use this fit to make precise predictions of 
other moments as a crosscheck. Bauer and Trott r i40l ex- 
amined fractional moments of the lepton spectrum with a 
variety of cuts, and found certain moments that were very 
insensitive to nonperturbative effects. These moments may 
therefore be predicted with small uncertainties using the 
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values of A and Ai determined from other moments, pro- 
viding a stringent test of the HQE for inclusive decays. 
Later, these moments were measured by CLEO[ .35J . and 
were found to agree beautifully with the theoretical predic- 
tions: 



D4 = 



.6GeV 



il.SGeV ^ 



172.3 rfr 

l.bGeY e dE, 



p2.9 cir 
SGeV { dEi 



0.5190: 
5193: 



0.0007 
0.0008 



0.6034 + 0.0008 
6036 + 0.0006 



(T) 
(E) 

(T) 
(E) 



(18) 



(where "T" and "E" denote theory and experiment, respec- 
tively). 

It should be noted, however, that there is currently poor 
agreement between the HQE and the first hadronic invari- 
ant mass moment, measured by BABAR with different lep- 
ton energy cuts, as shown in Fig. |3 This is a worrisome 
feature which should be better understood. One possi- 
ble resolution is to note that the BABAR measurement is 
not completely model independent, but depends on the as- 
sumed spectrum of excited D resonances. In the model 
used, there is no contribution to the B semileptonic width 
from excited D states with masses below ~ 2.4 GeV. This 
is an assumption which is in conflict with the results of the 
HQE, as noted a number of years ago by Gremm and Ka- 
pustin r l?n . 



may be obtained r i32l by assuming that all rate which does 
not go to the D or D* goes to the D** . This gives the in- 
equality 



{sH-mi,) = 0.051— +0.23 — + ... (19) 
\ n Mb I 

> TD"{m]y„ - m\,) + To-imo. - m\,) + Toim^, - m^) 



where Fx = F(B — » X£v)ITs.l. is the fraction of the 
semileptonic width which goes to the final state X, and the 
first line gives the first two terms in the HQE. Combining 
the HQE prediction with the measured ratio of the semilep- 
tonic B ^ D anAB D* widths r .42J 



0.31(1 -Fo-O, To. = 0.69(1 -Fz,..) 



(20) 



and taking mo" - 2.45 GeV leads to the upper bound 

Td" < 0.22 (21) 

which is in conflict with the experimental semileptonic 
branching fraction to excited states of ~ 0.35. Thus, if 
the prediction for the first moment of the hadronic invari- 
ant mass spectrum is valid, a non-negligible fraction of 
the semileptonic B width must be to excited D states with 
mx < 2.45 GeV, which would also bring the BABAR re- 
sults in better agreement with theory 1 1431 . It will be inter- 
esting to see how this situation evolves. 

3 Conclusions 



1 

0.8 
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0.4 
0.2 
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^ ^ I I 
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Figure 4. Comparison of the BABAR measurements of the 
hadron invariant mass spectrum vs. the lepton energy cut (black 
squares), and the HQE prediction not including BABAR hadronic 
mass data (red triangles), from Ref. [ iSSJ . 



Denoting the lowest excited state as D**, a lower bound on 
the first moment of the hadronic invariant mass spectrum 



The heavy quark expansion (HQE) has proven very suc- 
cessful at giving a model-independent description of inclu- 
sive B decays. Theory and experiment are now at the stage 
that a precision determination of \Vub\ is possible from in- 
clusive decays. The challenge is to incorporate kinematic 
cuts that exclude b —> c decays without introducing large 
uncertainties (theoretical or experimental). Cutting on the 
lepton invariant mass or an optimized combination of 
and the hadronic invariant mass mx gives a result that 
is insensitive to the nonperturbative light-cone distribution 
function fik+), at the expense of a difficult experimental 
measurement. Cutting on mx or the energy of the charged 
lepton Ee is easier experimentally, but introduces depen- 
dence on the nonperturbative parton distribution function 
f{k+). At leading order in 1 jniB, f{k+) may be determined 
from B —> Xsj decays, but there are potentially large sub- 
leading corrections to this relation (at least for the cut on 
Et) which may limit the ultimate precision of this method. 

Spectral moments from semileptonic b ^ c and radiative 
b s decays are now being used to study the HQE at the 
0(1 /m^) level, and to determine the values of the hadronic 
matrix elements which arise in the expansion. Global fits to 
a variety of moments now allow Vd, to be determined with 
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an uncertainty at the 2% level, and a short-distance nib with 
an uncertainty at the 100 MeV level. These determinations 
are not yet limited by theory, and so these determinations 
are likely to improve in the future. 



References 

1. F. Bartelt et al, CLEO Collaboration, Phys. Rev. Lett. 
71 (1993)4111; H. Albrecht er a/., Argus 
Collaboration, Phys. Lett. B255 (1991) 297. 

2. B.H. Behrens et al, CLEO Collaboration, 
iiep-ex/9905056 

3. P Ball and R. Zwicky, JHEP 01 10 (2001) 019. 

4. A. Kronfeld,hep-ph/0010074 

5. J. Chay, H. Georgi, and B. Grinstein, Phys. Lett. B247 
(1990) 399; M. Voloshin and M. Shifman, Sov. J. Nucl. 
Phys. 41 (1985) 120. 

6. I.I. Bigi et al, Phys. Lett. B293 (1992) 430; Phys. Lett 
B297 (1993) 477 (E); I.I. Bigi et al, Phys. Rev. Lett. 
71 (1993)496. 

7. A.V. Manohar and M.B. Wise, Phys. Rev. D49 (1994) 
1310. 

8. B. Blok et al, Phys. Rev. D49 (1994) 3356. 

9. A.H. Hoang, Z. Ligeti, and A.V. Manohar, Phys. Rev. 
Lett. 82 (1999) 277; Phys. Rev. D59 (1999) 074017. 

10. For a review, see A. El-Khadra and M. Luke, Ann. 
Rev. Nucl. Part. Sci. 52 (2002) 201. 

11. N. Uraltsev, Int. J. Mod. Phys. A14 (1999) 4641. 

12. V. Barger, C. S. Kim and R. J. Phillips, B Decays," 
Phys. Lett. B251,(1990) 629. 

13. A.F Falk, Z. Ligeti, and M.B. Wise, Phys. Lett. B406 

(1997) 225; I. Bigi, R.D. Dikeman, and N. Uraltsev, 
Eur. Phys. J. C4 (1998) 453. 

14. R. D. Dikeman and N. Uraltsev, Nucl. Phys. B509 

(1998) 378. 

15. P Abreu et al [DELPHI CoUaboration], Phys. Lett. 
B478 (2000) 14. 

16. M. Neubert, Phys. Rev. D49 (1994) 3392; D49 (1994) 
4623; T. Mannel and M. Neubert, Phys. Rev. D50 
(1994) 2037; I.I. Bigi et al. Int. J. Mod. Phys. A9 
(1994) 2467. 

17. F De Fazio and M. Neubert, JHEP06 (1999) 017. 

18. C.W. Bauer, Z. Ligeti, and M. Luke, Phys. Lett. B479 
(2000) 395. 

19. G. Buchalla and G. Isidori, Nucl. Phys. B525 (1998) 
333. 

20. M. Neubert, JHEP 0007 (2000) 022. 

21. C. W. Bauer, Z. Ligeti and M. E. Luke, Phys. Rev. 
D64 (2001) 113004. 

22. 1. 1. Bigi and N. G. Uraltsev, Nucl. Phys. B423 (1994) 
33; 

M. B. Voloshin, from inclusive semileptonic B 
decays," Phys. Lett. B515 (2001) 74. 

23. A. F. Falk, E. Jenkins, A. V. Manohar and 

M. B. Wise, decays," Phys. Rev. D49 (1994) 4553. 



24. A.K. Leibovich, I. Low, and I.Z. Rothstein, Phys. 
Rev. D61 (2000) 053006; Phys. Lett. B486 (2000) 86. 

25. M. Neubert, Phys. Lett. B513 (2001) 88. 

26. A. Bornheim et al [CLEO Collaboration], Phys. Rev. 
Lett. 88(2002) 231803. 

27. C. W. Bauer, M. E. Luke and T. Mannel, l/m(b)," 
|3-Xiv:hep-ph/0 102089 

28. C. W. Bauer, M. Luke and T. Mannel, |V(ub)|," Phys. 
Lett. B543 (2002) 261. 

29. A. K. Leibovich, Z. Ligeti and M. B. Wise, Phys. 
Lett. B539 (2002) 242. 

30. M. Neubert, Phys. Lett. B 543 (2002) 269. 
farXiv:hep-ph/0207002 1 . 

31. M.B. Voloshin, Phys. Rev. D51 (1995) 4934; M. 
Gremm et al, Phys. Rev. Lett. 77 (1996) 20. 

32. A. F Falk, M. E. Luke and M. J. Savage, Phys. Rev. D 
53, 2491 (1996); Decays," Phys. Rev. D 53, 6316 

(1996) ; A. F Falk and M. E. Luke, cut," Phys. Rev. D 
57, 424(1998). 

33. A. Kapustin and Z. Ligeti, Phys. Lett. B355 (1995) 
318; C. Bauer, Phys. Rev. D57 (1998) 5611; Erratum 
ibid. D60 (1999) 099907; Z. Ligeti, M. Luke, A.V. 
Manohar, and M.B. Wise, Phys. Rev. D60 (1999) 
034019; A.L. Kagan and M. Neubert, Eur. Phys. J. C7 
(1999)5. 

34. D. Cronin-Hennessy ef aZ. [CLEO Collaboration], 
Phys. Rev. Lett. 87, 251808 (2001). 

35. R. A. Briere et al [CLEO Collaboration], 
determination of the heavy quark expansion 
parameters," arXiv:hep-ex/0209024 

36. B. Aubert et al (BABAR Collaboration), 
"hep^x/0207084 

37. DELPHI Collaboration, Contributed paper for ICHEP 
2002, 2002-070-CONF-605; 2002-07 l-CONF-604. 

38. C. W. Bauer, Z. Ligeti, M. Luke and A. V. Manohar, 
m(b)," Phys. Rev. D 67, 054012 (2003). 

39. M. BattagUa et al, B decays," Phys. Lett. B 556, 41 
(2003). 

40. C. W. Bauer and M. Trott, Phys. Rev. D 67, 014021 
(2003). 

41. M. Gremm and A. Kapustin, Phys. Rev. D 55, 6924 

(1997) . 

42. K. Hagiwara et al [Particle Data Group 
Collaboration], Phys. Rev. D 66, 010001 (2002). 

43. see V. Liith, "Measurements of — > Xc(v Decays", 
talk given at ICHEP 2002. 



